Extremal statistics in the energetics of domain walls 
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We study at T = the minimum energy of a domain wall and its gap to the first excited 
state concentrating on two-dimensional random-bond Ising magnets. The average gap scales as 
l^Ei ~ L^fiN,), where f{y) ~ [lny]-l''^ 9 is the energy fluctuation exponent, L length scale, and 
the number of energy valleys. The logarithmic scaling is due to extremal statistics, which is 
illustrated by mapping the problem into the Kardar-Parisi-Zhang roughening process. It follows 
that the susceptibility of domain walls has also a logarithmic dependence on system size. 
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The energy landscapes of random systems are often 
assumed to be described at low temperatures by scaling 
exponents that follow from the behavior of the ground 
states. In the renormalization group (RG) language this 
means that temperature is an irrelevant variable. In most 
quenched random systems, the energy landscape contains 
many low-lying metastable minima separated by high 
barriers. Examples can be found in the realm of ran- 
dom magnets, the most famous one being of course spin 
glasses The dynamical behavior at finite tempera- 
tures, as a result of a temperature change or the appli- 
cation of an external field, will naturally depend on the 
associated barriers and energy differences between the 
minima. 

It is often assumed that energy differences or barri- 
ers between configurations (SE) relate to the length I 
involved by a scaling relation SE , where is an en- 
ergy fluctuation exponent. It measures the dependence 
of the first non-analytic correction to the ground state 
or free energy on the length scale. Here we show that 
for extended manifolds, or Ising magnet domain walls 
(DW) [equivalent to directed polymers (DP) in (IH-I)- 
dimensions] the energy difference between the ground 
state energy and the next state ( "first excited state" ) fol- 
lows from extremal statistics. This is due to the fact that, 
usually, one can assume that the energy landscape, at 
large enough scales, consists of many independent valleys. 
Finding the gap between the minimum and the second- 
most favorable state is then a straightforward extremal 
statistics problem as is the simpler one of the minimum of 
all the independent valley energies. The extreme statis- 
tics leads to logarithmic factors in the gap and minimum 
energies, which we show also by numerical calculations. 
The same result can be applied also to other disordered 
systems, where the energy landscape of DW can be re- 
duced to a one-dimensional form. We also interpret the 
results in the language of kinetic roughening, since DP 
maps into the Kardar-Parisi-Zhang (KPZ) equation of 
growth . Finally, as an application we show that the 
extremal statistics scaling shows up in the susceptibility 
of DW. 

Here we consider elastic manifolds at T = with 



quenched short-range, e.g. point-like defects, random- 
ness and in 0? = [D ~\-n),n = 1 dimensions, in which D is 
the dimension of the manifolds and d is the dimension of 
their embedding space. The continuum Hamiltonian for 
such an elastic manifold is 



n = 



-{Vz{^)Y + Vr{^,z) 



d^ 



X, 



(1) 



where z{'x.) is the height of the interface and x is the 
D dimensional internal coordinate of the manifold. The 
first term in the integrand is the clastic contribution 
with the corresponding surface stiffness F of the interface 
and the second term comes from the random potential. 
For random manifolds we use quenched random bond 
(RB) disorder, which means that the random poten- 
tial is delta-point correlated, i.e., (V^(x, z)V^(x', z')) = 
2I?(5(x — x')(5(z — z). The geometric behavior of the man- 



ifold is characterized by uP' 



z(x) 



where L is the linear size of the system and C is the cor- 
responding roughness exponent. At low temperatures in 
(1 + 1) dimensions, due to the equivalence of DP in ran- 
dom media |^,|| to the KPZ equation, the exact rough- 
ness exponent reads Q = 2/3 1^-^. In higher dimen- 
sions functional RG gives the approximate expression 
C ~ 0.208(4 - D) |] for RB DW. Since the width of 
a manifold grows as L*" it is expected that the number of 
independent valleys is proportional to ~ L^/L'^. 
At T = the total average minimum energy {Eq) of an 
elastic manifold equals its free energy and grows linearly 
with the manifold area and its fluctuations scale as 
A£; = ((£^0 - (£^0))^)^ - i^ where ^2C + D-2^. 

Let us now derive analytically the scaling of the "ex- 
treme statistics" contributions to the lowest minimum, 
Eq, and the gap between two lowest minima, Ai?i = 
El — Eq. We consider the case of many independent val- 
leys in the landscape > 1, which means that the DP 
can have an arbitrary starting or end point, and that 
Lz > L^. For the "single valley" boundary condition 
case (one end of the manifold fixed) , it is known numeri- 
cally that near its mean the distribution is Gaussian 
Hence we draw the energies E from the distribution 



1 



P{E) ^kcxp<- 



\E-{E)\ 
AE 



(2) 



where (E) ^ is the average energy of the manifold 
and AE ~ measures its fluctuations and k normalizes 
the integral so fc The exponent is not con- 

stant 1^,^. Near the peak rj — 2. In the low energy tail 
numerical simulations indicate that 77_ « 1.6, while in 
the high energy tail the best estimate is 77+ « 2.4 At 
this stage we allow rj to be variable but note that it is the 
behavior near the mean and the low energy tail which is 
the most important in this calculation. In a system with 
Nz ^ Lz/L'^ independent local minima the probability 
that the global minimum has energy E is given by, 



Lr,AE)^NzP{E){l-Ci{E)} 



(3) 



where Ci{E) = j_^P{e)de [0. The gap AEi follows 
similarly. Its distribution, Gjv^ {AEi , E) is given by 



Gn^ iAEi,E) = NziNz - l)P{E)P{E + AE^) 
{l-Ci(£; + ASi)}^--2. 



(4) 



Gm^ {AEi, E) is the probability that if the lowest energy 
manifold has an energy i?, then the gap to the next low- 
est energy level is AEi . The average value of the global 
minimum is given by 



{Eo) 



ELNAE)dE, 



(5) 



which is not analytically integrable. The typical value 
of the lowest energy may be estimated using an extreme 
scaling estimate. It follows from the fact the term the in- 
side { } in (||) becomes unity if Ci is small enough. It has 
proven useful in other contexts, for example breakdown 
of random networks, and reads here [Tl| ], 



l/kNzP{{Eo)) « 1 



which yields. 



{Eo) « (E) - AEiHNz)}'/^, 



(6) 



(7) 



where AE L^. 

To estimate the typical value of the gap, we use simi- 
larly to (|) 

l/k''Nz{Nz - l)Pi{Eo))Pi{Eo) + (ASi)) « 1, (8) 
which with and the fact that \{AEi)\ <C \{Eo) \ yields, 
AE'J AE 



{AE,) 



rj{{E) - (i?o))"-i ^ 77{ln(iV,)}("-i)/"- 



(9) 



We thus find that in addition to the usual sample to sam- 
ple variations in the energy {AE ~ L^) there is a slow 
reduction in the gap which scales as {ln(A^2)}~^''^^)/'', 
provided A^^ > 1. Our case is closely related to the weakly 



broken replica symmetry of DP, see also Ref. [ [L3[ , 
where the relation between replica methods and extremal 
statistics is discussed. 

The (1-1-1) dimensional DW maps, in the continuum 
limit, to the KPZ equation by associating the minimum 
energy of a DW with the minimum arrival time ti = Eq 
of a KPZ-surface to height h. The connection is illus- 
trated in Fig. |l| in the limit of many valleys > 1. 
The minimal path of the DW with the endpoint z{L) 
equals the path by which the interface reaches h = L, 
at location xi = z and at time ti = Eq. Thus ti 
attains a logarithmic correction, from Eq. of size 
-hl^{\n{Lz/h^/'')y/'^, where (3 = 1/3 and z = 3/2 are 
now the roughening exponent and dynamical exponent 
of the KPZ universality class Consider now the sec- 
ond smallest arrival time t2- In the KPZ language of 
DP, if the path X2{t') that gives ^2 is completely inde- 
pendent of the Xi{t') that results in ti then t2 and X2 
are related to a separate, independent valley of the DP 
landscape. The difference At = t2 — <i then equals AEi 
of DW, and likewise obeys extremal statistics, so that 
At - ft.^[ln(Ljft,i/^)]-(''-i)/''. For growing surfaces this 
limit is the early stages of growth, in which the correla- 
tion length ^ <C L2, and therefore the arrival times, or 
DW energies, are independent. 

In order to check the scaling behavior of the gap en- 
ergy (^, we have done extensive exact ground state cal- 
culations of elastic manifold in the two dimensional (2D) 
spin-half RB Ising model, i.e. we take a nearest neigh- 
bor Ising model with random but ferromagnetic couplings 
Jij > 0. Calculations are performed by varying both the 
parallel length L and the height of systems oriented 
in the {10} direction. The DW is imposed by antiperi- 
odic boundary conditions in the z-direction at z = and 
z — Lz- The elastic manifold is the interface, which 
divides the system in two parts one containing up spins 
and the other containing down spins. At T = the prob- 
lem of finding the ground state DW is a global optimiza- 
tion problem which is solved exactly using a mapping to 
the minimum-cut maximum-flow problem. The so-called 
push-and-relabel method solves this problem efhciently 
and has been extensively discussed elsewhere p"^-p^. 

In order to control the average number of the minima 
(A'^) ^ Lz/L'' in a chosen system size, we set the ini- 
tial position of the interface zq in a fixed size window 
at height 'zq/Lz — const. If the ground state interface is 
originally outside the window, with room only for a single 
valley, it is neglected and a new configuration is created. 
After the original ground state is found, with its energy 
Eq , the lattice is reduced so that bonds in and above the 
window are neglected and the new ground state, its Ei 
and the corresponding gap energy AEi, is found. We 
studied at least A^ = 500 realizations of system sizes up 
to L = 300 and Lz = 500. Fig. || starts the discussion 
of the numerical data by showing how the ground state 
energy (Eq) behaves as a function of L and Lz- The 



2 



scaling result (|^) expects that the correction to the en- 
ergy follows a logarithmic dependence on N^, which is 
confirmed in the figure. Note that the extraction of this 
correction from the data requires an educated guess of 
how (E), the single valley energy, behaves with L. We 
have used an Ansatz (E) aL + b with the values of a 
and b demonstrated in Fig. ^ so that the exponent value 
ry = 2 corresponds to a Gaussian distribution. Due to the 
nature of the procedure it would probably be possible to 
obtain a reasonable fit for e.g. 77 = ry_ as well. 

For small sample sizes, < L'' the value of the en- 
ergy Eq is affected by confinement. Similarly, the gap is 
controlled by confinement effects in this limit. When 
is large there are many independent valleys and extreme 
statistics effects are important, hence we expect, 



(A^i (L,L,)) 



Lz < L, 
> L 

(10) 

where we have used Eq. (||) and ^ L^/L''. We at- 
tempt to collapse the data by using the reduced variables 
{AEi{L, Lz))JL^ versus L^/L'' for various L and L^. As 
seen in Fig. we find a nice agreement with the extreme 
scaling form, with the ratio (77— l)/ri = 1/2, i.e. by using 
a Gaussian distribution. 

Next we consider the relation of the extremal statis- 
tics to the susceptibility of these manifolds. In the D- 
dimensional case the susceptibility is defined by, 



/ dm 
h^o+ \ ah 



(11) 



where the change in the magnetization of the whole d 
dimensional system is calculated in the limit of the van- 
ishing external field from the positive side ||l^,|l^ and 
the brackets imply a disorder-average. We have recently 
shown that the general behavior follows from a level- 
crossing phenomenon, which involves an extra potential 
Vh{z) = hz, dependent on the height of the interface, in 
the Hamiltonian (|l|), and h is an applied external field 
to the manifold. In any particular configuration when h 
is varied the manifold position changes in macroscopic 
'jumps' iQ, the first one occurring at hi. 

One may write the susceptibility, Eq. (pi]), with the 
help of the probability distribution of the fields hi P{hi) 
in the form 



X = lim ( — — 

h^o+ \ Ah 



Azi 



lim P(hi), 
/1-.0+ 



(12) 



because the magnetization of a system m(/i) ~ z{h)/Lz, 
and since the distance in the jump between the min- 
ima (Azi) ^ Lz ||l^, independently of the sample- 
dependent hi. It is expected that a scaling form P{hi) ~ 
{hi) P {hi/ {hi)) applies, and that P remains finite in the 
limit hi — > 0. Next we compare the average suscepti- 
bility as a function of the number of valleys Nz to the 



conjecture that in the presence of the field the average 
gap for the original and excited state follows an extremal 
statistics form similar to Eq. (^). 

The simulations are done again using a fixed height 
window in which the original ground state without the 
field is found. After that the external field h is slowly ap- 
plied by increasing the coupling constant values J±{z) = 
J random + hz, where J_L is perpendicular to the z- 
direction, until the first jump is observed with the cor- 
responding hi and Azi. In order to find the scaling re- 
lation for the first jump field /ii, we make the Ansatz 
{AEi) — {hi)LLz, since the field contributes to the poly- 
mer energy proportional to {D = 1) and Lz ^ {Azi) 
is the difference in the field contributions hz to the en- 
ergy at finite h at different average valley heights zq, zi. 
Hence 



{hi{L,Lz))LLz^ L'f 



(13) 



where the scaling function /(y) — [ln{Lz/L'')]^'^^^''/^ . 
Fig. ^ shows the scaling function (13) with a collapse 
of {hi{L, Lz))L^~^ Lz versus Lz/L^ for various L and 
Lz which is again in good agreement the logarith- 
mic extreme scaling correction. Generalizing to arbi- 
trary dimensions one has the behavior of {hi{L, Lz)) ~ 
i^-^_L^ i[ln(LjL^)]-(''-i)/''. For the susceptibility, 
Eq. (13), one obtains, using {hi) for the normalization 
factor at P{hi — 0), 



X-L°-'Lz[HLz/Li)]<^^-'y^, 



(14) 



and in the isotropic limit, L cc Lz, the total susceptibility 
Xtot — L'^X becomes (when rj = 2) 



Xtot-L'^'+'-'ai-OHL)]'/'. 



(15) 



Notice that for most random manifolds 1 — C > with 
the exception of 2D random field Ising DW for which 
C ~ 1 at large scales |18 , and thus the susceptibility 
does not diverge |jl^ as the premise A^z > 1 does not 
hold in that case. If the condition A^^ > 1 is violated 
the extreme statistics correction disappears. Thus the 
extremal statistics of energy landscapes leads to a loga- 
rithmic multiplier in the susceptibility, Eq. ([T^), of the 
DW. This result differs from algebraic forms of scaling 
]l§, see also 

To conclude, we have considered the average energy 
differences or "gaps" in the energy landscape of (two- 
dimensional) elastic manifolds. An extremal statistics 
argument in the system geometry with many indepen- 
dent valleys shows that the ground state energy and the 
gap have logarithmic scaling functions, also reproduced 
with numerical studies. An illuminating connection can 
be made to Kardar-Parisi-Zhang noncquilibrium surface 
growth. Finally, we demonstrate that the gap scaling 
shows up in the susceptibility of random manifolds. This 



3 



might have imphcations for flux hne lattices in high- 
temperature superconductors, where a similar problem 
related to barriers has been analyzed with the aid of ex- 
tremal statistics [ElJ. 
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FIG. 1. The relation between DP and growing interfaces. 
KPZ interface is growing so that h increases and DPs in in- 
dependent valleys equal the n-th fastest arrival times of the 
interface to a prefixed height h, at Xn, at times t{xn) in a 
system with width Lz- The solid line describes the fastest 
polymer, which ends at x\. The dashed lines describe the 
next fastest polymers. 
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FIG. 2. The scaling of the ground state energy Eo 
as a function of scaled transverse system size L^/L'' for 
the system sizes L =100, 200, and 300. The line 
-0.41 + 0.53[ln(2.78L^/L^)]^/^ is a guide to the eye. We have 
subtracted the expected dependence of (E) from (Eq) (see 
the text). We use in Figs. §to ^ RB disorder with Jij.z £ [0- 
— 1] uniform distribution and Jij,x ~ 0.5. The number of 
realizations ranges from A'^ — 500 for L = 300, = 500 to 
N = 2000 for L = 200, = 600. 
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FIG. 3. The scaling function f{y) of the scaled disor- 
der-average of the energy difference {AEi)/L^ as a function 
of scaled transverse system size Lz/L'' for the system sizes 
L =100, 200, and 300, each with zo/L^ ~ const. 9 = 1/3, 
C = 2/3. The line has a shape f{y) = 0.23 ln(i/)"^/^. The 
configurations are the same as in Fig. B. 
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FIG. 4. The scaling function f{y) of the scaled disor- 
der-average of the jump field {hi)L^~^ as a function of 
scaled transverse system size L^/L^ for the system sizes 
L =100, 150, 200, 250 and 300, each with zo/L^ ~ const, 
e = 1/3, C = 2/3. The line has a shape f{y) = 0.41 \n{y)-^''^. 
Here the number of realizations ranges from = 500 for 
L = 300, Lz = 500 to Af = 2600 for L = 200, = 600. 



5 



